Defect generation in a spin-1/2 transverse XY chain under repeated quenching of the 

transverse field 
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We study the quenching dynamics of a one-dimensional spin-1/2 XY model in a transverse field 
when the transverse field h(— t/r) is quenched repeatedly between — oo and +oo. A single passage 
from h — > — oo to h — > +oo or the other way around is referred to as a half-period of quenching. 
For an even number of half-periods, the transverse field is brought back to the initial value of — oo; 
in the case of an odd number of half-periods, the dynamics is stopped at h — > +oo. The density of 
defects produced due to the non-adiabatic transitions is calculated by mapping the many-particle 
system to an equivalent Landau-Zener problem and is generally found to vary as 1/\/t for large r; 
however, the magnitude is found to depend on the number of half-periods of quenching. For two 
successive half-periods, the defect density is found to decrease in comparison to a single half-period, 
suggesting the existence of a corrective mechanism in the reverse path. A similar behavior of the 
density of defects and the local entropy is observed for repeated quenching. The defect density 
decays as 1/\/t f° r large r for any number of half-periods, and shows a increase in kink density for 
small r for an even number; the entropy shows qualitatively the same behavior for any number of 
half-periods. The probability of non-adiabatic transitions and the local entropy saturate to 1/2 and 
In 2, respectively, for a large number of repeated quenching. 

PACS numbers: 73.43.Nq, 05.70.Jk, 75.10.Jm 



I. INTRODUCTION 

Zero-temperature quantum phase transitions 1 ' 2 driven 
by quantum fluctuations have been studied extensively in 
recent years. In a quantum system, statics and dynam- 
ics are intermingled and hence a quantum critical point 
is associated with a diverging length scale as well as a 
diverging time scale called the relaxation time. The re- 
laxation time of a quantum system is the inverse of the 
minimum energy gap which goes to zero at the quantum 
critical point. In the proximity of a quantum critical 
point, the spatial correlation length £ grows as £ ~ \8\~ v 
and the characteristic time scale or the relaxation time 
£ T scales with £ as £ T ~ £ z , where S is the measure of 
the deviation from the critical point, and v and z are the 
critical exponents characterizing the universality class of 
the quantum phase transition. 

When a quantum system is swept across a zero- 
temperature critical point 1,2 by slowly varying a param- 
eter of the Hamiltonian at a uniform rate, the result- 
ing dynamics fails to be adiabatic, however slow the 
time variation may be. This is because of the diverg- 
ing relaxation time discussed above. There exist non- 
adiabatic transitions which eventually lead to defects in 
the final state. According to the Kibble-Zurek (KZ) 
argument^, the non-adiabatic effects dominate close to 
the critical point where the rate of change of Hamilto- 
nian is of the order of the relaxation time of the quantum 
system. The KZ analysis predicts that when a parame- 
ter of the Hamiltonian is quenched at a uniform rate as 
t/r through the critical point, the density of defects in 
the final state shows a power-law behavior with the time 
scale of quenching r given by l/T ud ^ uz+1 \ Following 
recent experimental studies^ of non-equilibrium dynam- 



ics of strongly correlated quantum systems, there is an 
upsurge in the theoretical investigation of related non- 
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and random models*., and models with a 



gapless phased. A generalized scaling relation for the 
defect density in a non-linear quench across a quantum 
critical point has also been proposed^. 

The quenching dynamics of a spin-1/2 XY chain in a 
transverse field, when either the transverse fielc£ or the 
interaction^ is quenched from — oo to +oo at a uniform 
rate t/r, has been studied extensively in recent years, and 
the defect density is found to obey the KZ prediction. 
It should also be noted that recent studies of the two- 
dimensional Kitaev model indicate a generalization of the 
KZ prediction when quenched along a critical surface 12 . 

In this work, we investigate a situation where the trans- 
verse field is quenched back and forth across the quantum 
critical points from — oo to +oo and again from +oo to 
—oo and so on, with the functional form of the transverse 
field being given by h — ±t/r. According to the Kibble- 
Zurek argument, the system fails to evolve adiabatically 
near the quantum critical points, where the relaxation 
time is very large, which results in the production of kinks 
(in this case, oppositely oriented spins in the final state). 
In our notation, nil) corresponds to the defect density 
in the final state after I passages through both the Ising 
quantum critical points to be defined below. Thus even 
values of I signify I half-periods and the transverse field is 
brought back to the initial value h — > — oo; for odd values 
of I, the field in the final state tends to +oo. Clearly, the 
case of I = 1 has been extensively studied earlieri&. In 
all the cases, the initial value of h is large and negative 
so that all the spins are down, i.e., oriented antiparallcl 
to the z-axis. If the dynamics were adiabatic for the en- 
tire span of time, one would expect all spins to be down 
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(up) in the final state following an even (odd) number 
of half-periods of quenching, and the defect is defined 
accordingly. 



II. THE QUENCHING SCHEME AND RESULTS 

The Hamiltonian of our model is given by 
H = - ^ ( J *«+i + J v«+i + h<), (!) 

n 

where the cr's are Pauli spin matrices satisfying the usual 
commutation relations. The strength of the transverse 
field is denoted by h, and J x — J y is the measure of 
anisotropy; J x , J y and h are all non-random variables. 
The Hamiltonian in Eq. (fTJ) can be exactly diago- 
nalized using the Jordan- Wigner transformation which 
maps a system of spin-l/2's to a system of spinless 
fermion o 13 ' i 15 . Diagonalizing the equivalent fermionic 
Hamiltonian in terms of Bogoliubov fermions, we arrive 
at an expression for the gap in the excitation spectrum 
given b y 13 i 15 



final time is given byA£ p k = |Cifc(+oo)| 



2lr7 , where 



[h 2 + J 2 + J 2 +2h(J x + J y ) cos k+2J x J y cos2/c] 1/2 .(2) 

The gap given in Eq. |2]) vanishes at h = ^{J x + Jy) 
for wave vectors k = and tt respectively, signaling a 
quantum phase transition from a ferromagnetically or- 
dered phase to a quantum paramagnetic phase known as 
the "Ising" transition. On the other hand, the vanishing 
of the gap at J x — J y for \h\ < (J x + J y ) at an order- 
ing wave vector ko = cos (h/2J x ), signifies a quantum 
phase transition which belongs to a different universality 
class from the Ising transitions between two ferromag- 
netically ordered phases. In our quenching scheme, the 
system will be swept across the Ising critical points only. 

Let us first briefly recall the case with I = 1— When 
projected to the two-dimensional subspace spanned by 
the state vectors |0) (empty state) and \k,— k) (two 
fermion state), the Hamiltonian takes the form 



H k (t) 



- [h + {J a 
h+{J 2 



Jy) COS k] I 2 

-J y )cosk i(J x — J y ) sink 
J y ) sin k —h — (J x + J y ) cos k 



where I2 denotes the 2x2 identity matrix. Therefore, the 
many-body problem is effectively reduced to the prob- 
lem of a two-level system, with the two levels being the 
states |0) and \k, — k). The general state vector tpkif) 
evolving in accordance with the Schrodinger equation 
idtipk(t) — Hk(t) ipk(t) can be represented as a linear 
superposition \ipk(t)) = Cifc(i)|Q) + C 2k (t)\k, -k), with 
the initial condition C\k{— 00) = 1 and C 2k {— 00) = 0. 
The off-diagonal term of the projected Hamiltonian, A = 
( J x — J y ) sin k, represents the interaction between the two 
time-dependent levels £1,2 = ±[h(t) + (J x + J y )cosk]. 
The Schrodinger equation given above is identical to the 
Landau-Zener (LZ) problem of a two-level system^ and 
therefore, the transition probability of excitations at the 



7 = A 2 /a, so that in the present casepk(l = 1) = Pk(l) — 
e --KT(Jv-Jy) sin yhis immediately leads to an expres- 
sion for the density of kinks or down spins n in the final 
state 



<-x <ik 

n(l = 1) = n(l) = / — p k 



1 



T^\fr \ Jx - Jy\ 



(3) 



in the limit of large r. The 1 / yfr decay of the defect den- 
sity is in accordance with the Kibble-Zurek prediction^. 

We shall now generalize the quenching dynamics to the 
case 1 = 2, when the system is brought back to the state 
with h — > —00 from the final state of the case 1 = 1, 
with the initial condition given by |Cifc(oo)| 2 = e~ 2 '* 1 
and |C2fe(oo)| 2 = 1 — e~ 27r7 . With a view to estimating 
the non-adiabatic transition probability, we consider two 
time-dependent states |1) and |2) with energy ei(t) and 
£2(^)7 respectively where £i(t) — ^(t) = at, a being a 
constant and the time-independent coupling between the 
states being A. Let us also assume that the time t goes 
from —00 to +00 (the forward path in Fig. 1). 



Forward Path 

> 


^ |l(oc)> 
2A 




Time 

V Reverse Path 




■■■■■'--....„ ^ < 

2^> 



FIG. 1: The time-dependent energy levels of the Landau- 
Zener Hamiltonian. The minimum gap is 2A. 

Defining a general state vector as \^f(t)) = Ci(z:)|l(<)) + 
C 2 (t)\2(t)), where \d{t)\ 2 {\C 2 {t)\ 2 ) is the probability of 
the state |1)(|2)) at time t, we can recast the Schrodinger 
equation in the form 



d 2 - 



1 



1 ^U(z) + (r + ---z 2 )U(z)=0, 



(4) 



where z = e ^aSt, r = iA 2 /a, and U(t) = 
U(z) is related to C 2 (t) through the relation U(t) = 
C 2 (t)e l S t -<~ dt '^ {t ' ) ^ fL dt '(^')-^(t')) . Focusing on the 
special case with \C\{— 00) | = 1 and |Ca( — oo)| = 0, 
we have U(z) = AD- r -i(—iz) as the particular solu- 
tion of the Weber's differential Eq. (4) 17 i 18 . We re- 
mark that the axis in the z plane which corresponds to 
t, is along e^J 7 ™ for t > 0, and along ei m for t < 0. 
The constant A is determined from the initial condition 
|Ci(— 00) I = 1 by taking the asymptotic expansion of 
D- r -i(—iz) along e^" (or t — > —00 limit). This fi- 
nally yields D- r -x(-iRei ni ) w e -f («+i)» e -^ 2 /4 jR -«-i j 
where R = y/\a~\\t\ and we find \A\ = ^^e-^^. The 
solution in the limit R — > +00 (or t — > +00) along z = 
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leads to the final expression |Ci(oo)| 2 = e~ 27r7 
(a result already used for the case 1 = 1). The parameter 
7 = A 2 1 a depends on the magnitude of the slope a of 
the two approaching states and the interaction A. Using 
the above results, we therefore get a recursive relation 
for the probabilities after the (I + l)-th quenching (with 
I > 0) given as 



|Ci(-(-l) 



< +1 oo)| 2 



e - 2 ^|Ci(-(-l) z oo)| 2 
(l-e- 2 ^)|C 2 (-(-l)'oo)| 2 , (5) 



|C 2 (-(-l)' +1 co)| 2 = (l-e-^ldH-l^oo)! 2 



\C 2 (-(-l) l oo)f 



(G) 



Eqs. (5) and (6) look incomplete at first sight 
because they do not contain any cross terms like 
Ci(±oo)C2(±oo)* or C2(±oo)Ci(±oo)*. However, this is 
because at t — > ±oo, the coefficients C\(t) and C 2 (i) vary 
rapidly with time as exp[±(i/fi.) J* dt'E(t')^. Hence, the 
two cross terms given above vary rapidly with the initial 
time (which is going to +oo or — oo depending upon the 
number of repetitions as explained above) and indepen- 
dently of each other for different values of k. Their con- 
tribution to the defect density therefore vanishes upon 
integration over k due to the presence of terms like 
cos(Tcos/c), T being the time at which the terms are cal- 
culated (which is ±oo in our case). On the other hand, 
the terms given in Eqs. (5) and (6), namely |Ci(±oo)| 2 
and |C 2 (±oo)| 2 have no such rapid variations, since any 
arbitrary large phase in C\ (±oo) cancels the exactly op- 
posite phase appearing in Ci(±oo)* (with a similar argu- 
ment holding for C2(±oo). Thus Eqs. (5) and (6) follow 
from an exact solution of the Landau-Zener problem with 
general initial conditions Ci(— oo) and C 2 (— oo), once we 
use the fact that the phases of these initial amplitudes 
are rapidly varying and are therefore uncorrelated with 
each other. This is also the reason why we can consider 
the state obtained after one or more quenches as a mixed 
state with an entropy given by the expression in Eq. (9) 
below. 

Using the above results, one finds the probability of a 
non-adiabatic transition to be 



P*(2) = |C 2 (-oo)| 2 = 2e- 2 ^(l - e- 2 ^). 



(7) 



The transition probability pk (2) (Fig. 2) is maximum for 
k = tt/2 for small values of r, whereas for higher values 
of r, there are two maxima symmetrically placed around 
1 2 which gradually shift 



k — 7T / 2 clt /en — oil. i t — ~ . 

I U V ITTiJ^-Jy) 2 

to k = and 7r for very large r. In the limit of r — > +oo, 
there is no non-adiabatic transition as the dynamics is 
perfectly adiabatic. The observed behavior of Pfc(2) also 
corresponds to the existence of an inherent time scale tq 
(for fco = 7r/2) in the problem which separates the regions 
of small and large r. It is interesting to note that the 
non-adiabatic transition probability for the forward path 
attains the minimum value 7 of 1/2 at the wave vector 



k = n/2 for r = t = (ln2)/[ir(J x — J y ) 2 ], while for the 
reverse path, although Pfc 0=w /2,r (2) is once again 1/2, 
this is the maximum possible value of pk(2). 



- T = T„=0.22 



-> X = 600 



X= 1 



x=o.oi \ 




FIG. 2: Variation of p*,(2) with k for different values of r as 
obtained analytically for J m — J y = 1. 

The kink density n(2), i.e., the density of the up spins 
in the final state with h — > — oo, is again related to p/c(2) 
as n(2) = - dkpk(2) (Fig. 3 (a)). In the limit of large 
r, we find that 



n(2) 



TT-Jt(J X — Jy 



;0 



(8) 



Clearly, the magnitude of the defect density after a full 
period (I = 2) case, is smaller than the / = 1 case given 
in Eq. (3) in the limit of large r, and also in the limit of 
small r when the defect density is maximum for 1 = 1. 
This establishes the existence of a corrective mechanism 
during the reverse quenching process, arising from the 
fact that the maximum possible value of j?/c(2) is 1/2, 
which makes the area under the curve of pk vs k smaller 
for the / = 2 case as compared to the I = 1 case. The 
defect density n(2) (see Fig. 3) attains a maxima around 
r ~ 2tq; eventually there is a 1/y/r decay in the asymp- 
totic limit of t. We can explain the n vs t behavior in 
the following way: in the limit of small r, the system 
fails to evolve appreciably, and always remains close to 
its initial state for both I = 1 and 1 = 2, whereas for 
large r, the system evolves adiabatically at all times ex- 
cept near the quantum critical points. In either situation 
the non-adiabatic transition probability pk (2) and hence 
the density of defects n(2) are small though it is larger 
in the intermediate range of r. 

Although the final state after a full period is a pure 
state, locally it may be viewed as a mixed state described 
by a decohered reduced diagonal density matrix with 
elements^ pk and 1 — pk- The Von Neumann entropy 
density of the final state is given by 



* dk 



PkHVk) + (l-pfe)ln(l-pfe)]. (9) 



s(2) shows a similar behavior with r as n(2) (Fig. 3(b)), 
and the behavior can be justified along the same line of 
arguments as given above. 
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> n (4 ) 
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(see Fig. 3 (a)) also unravels other interesting aspects of 
the repeated quenching dynamics: (i) For even values of 
I, the characteristic scale t* at which the defect density is 
maximum shifts to higher values of t (from t* ~ 2to for 
I = 2). At the same time we find that n(l) > n(l — 2) for 
all values of r so that it eventually saturates to 1/2 for all 
r for asymptotic I. (ii) Similarly for odd I, n(l + 2) < n(l) 
for small r; for higher r, n(l + 2) exceeds n(l) following 
a crossover around t ~ r*. Once again n(l) saturates to 
1/2 for a large number of repeated quenches. The be- 
havior of the entropy density as a function of r for any 
general value of I is shown in Fig. 3 (b). It follows from 
Eq. (9) that for any given r, the local entropy density in- 
creases with I; in the limit I — > oo, Pk(l) tends to 1/2 and 
therefore the local entropy density tends to its maximum 
value of In 2. 



FIG. 3: (a) Variation of n(l), n(2), n(3) and n(4) with r, 
and (b) variation of s(l), s(2), and s(3) with r as obtained by 
numerically integrating Eqs. (9) and (10) with J x — J y — 1. 
In (a), n(4) denotes the defect density as obtained from the 
analytical expression in Eq. (11). In the limit of large r, the 
numerical results match perfectly with the analytical results. 



For an arbitrary number of repeated quenching, the 
non-adiabatic transition probability Pfc(Z) = (1 — e -2 *' 7 ) — 
(l-2e~ 27r7 )(l-p fe (Z-l)). We can simplify this to obtain 



Pk(l) 



I (1 - 2e- 2 ^) 1 



(10) 



Eq. (10) shows that pu{l) increases with I for even values 
of I, while for odd I, this variation depends on the value of 
k. However, in the asymptotic limit of I, pk(J) saturates 
to the value 1/2 for all k which implies that every spin 
is up or down with an equal probability in the final state 
following a large number of quenches. 

Using pk(l) given above for any I, we find that the 
defect density scales as \j\pr in the asymptotic limit of 
r and can be put in the form 



n(l) = 



1 



2tt{J x - Jy)^ S w \(l - w )\ ^ 



n 



Interestingly, using Eqs. (10) and (11), for two succes- 
sive values of I, we find that n(l + 1) < n(l) for small 
values of r if I is odd, a fact that once again emphasizes 
the corrective mechanism in the reverse path mentioned 
already. A similar result is also obtained for large r for 
smaller values of I, as shown in Fig. 3 (a). On the other 
hand, n(l + 1) is always greater than n(l) for even I. We 
also find that n(l + 1) — » n(l) in the asymptotic limit of 
I and r for both even and odd values of I. Moreover, for 
any odd I, the kink-density n(l) decreases monotonically 
with t, while for even I, n(l) attains a maxima around 
To; in the limit r — > oo, n(l) decays as 1/\/t in both 
cases. A close inspection of the variation of n(l) with r 



III. CONCLUSION 

In conclusion, we have studied the defect generation 
and entropy production in a transverse XY spin-1/2 
chain under repeated quenching of the transverse field 
between — oo and +oo . We have employed a generalized 
form of the Landau-Zener transition formula along with 
the concept of uncorrelated initial phases of the probabil- 
ity amplitudes so that the cross terms appearing in the 
recursion relation of probabilities vanish under integra- 
tion over momentum. Using the non-adiabatic transition 
formula thus obtained, we evaluate the defect density in 
the final state after an arbitrary number of quenches. 

Our results show that the defect density satisfies n(l) > 
n(l + 1) for small r, if I is odd; this points to the existence 
of a corrective mechanism in the reverse path. The re- 
sults obtained by numerical integration of the transition 
probability and by using the analytical expression given 
in Eq. (11) match perfectly in the limit of large t. The 
entropy density, however, is found to increase monoton- 
ically with the number of repetitions, showing that the 
local disorder of the system increases monotonically with 
I, irrespective of the behavior of the kink density. For an 
odd number of repetitions, we observe a monotonic de- 
crease of the kink density with t, as seen previously for 
the widely studied I — 1 case. For even /, on the other 
hand, n(l) grows for small r but eventually decreases as 
1/\/t in the large r limit, attaining a maxima at an in- 
termediate value of r = t*; t* shifts to higher values of r 
with increasing I. The difference in the behaviors of the 
defect density for an even and odd number of repetitions 
originates from the fact that the system is expected to 
come back to its initial ground state after an even value 
of I for a perfectly adiabatic dynamics. 

As mentioned above, the local entropy increases after 
each half-period of quenching. In the limit I — > oo, the 
entropy s(l) eventually saturates to its maximum possible 
value of In 2 for all t while the non-adiabatic transition 
probability Pk(l) approaches 1/2 for all k. This result 
suggests that a spin remains up or down with the same 
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probability, irrespective of the applied field after a large 
number of quenches. Using Eq. (10), one can also define 
a characteristic number of repetitions l*(k,r) so that for 
I > l*(k, r), the transition probability pk(l) ~ 1/2; it can 
be shown that l*(k, r) attains a minima at an interme- 
diate value of r. We conclude by noting that the defect 
production due to repeated quenching can be studied in 
systems of atoms trapped in optical lattices^, quantum 
magnets^ and spin-one Bose condensates^. 
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